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Compton Observatory Observations of Clusters of Galaxies and 
Extragalactic Radio Sources. 

Science Summary. 

This task involved the investigation of the emission of clusters of 
galaxies, particularly those which contain extended radio emission, in 
the gamma-ray region of the spectrum. Observations were made of several 
clusters using the Compton Observatory EGRET instrument. For each 
cluster a measured flux or upper limit on the gamma-ray flux was 
obtained. In only one case, Abell 2199, was there a significant measured 
flux. This source is spatially confused with a known blazar in the field 
of view. The observation is consistent with all emission being from the 
blazar . 

Essentially, this program has led to an improvement of an 

order of magnitude on the upper limits for the gamma-ray luminosity 

of clusters of galaxies. This result has important implications 

in the studies of the metagalactic cosmic ray background, the 

composition of the unseen matter in clusters of galaxies and 

the mechanisms by which some clusters of galaxies generate synchrotron 

radio emission over megaparsec scales. Preliminary results of this 

analysis has been published in the Proceedings of the Second 

Compton Symposium (AIP Conference Proceedings 304, p669, 1993) and in 

conjunction with later observations at the High Energy Astrophysics 

Division meeting of the AAS in November 1994. A final paper, funded 

separately, is being prepared for submission to the Astrophysical Journal . 

Observations and Analysis: 

Clusters observed specifically for the program included: 

Abell 426 
Abell 2162 
Abell 2199 
Abell 2255 
Abell 2256 
Abell 2319 

The following clusters were analyzed using Compton Observatory 
data : 

Abell 400 
Abell 496 
Abell 1060 
Abell 1367 
Abell 1656 
Abell 2022 
Abell 2052 
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These observations were processed using a combination of both instrument 
team software and software developed under this task. Analysis tasks 
included adding together data from overlapping viewing periods, 
developing tools to measure source and background counts, and 
determining confidence limits on detection. The EGRET team Like program 
was used in the final step of the analysis but all other software 
was developed under this program. 
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Abstract. 


Clusters of galaxies may emit high-energy gamma rays due to a number of emission 
mechanisms including cosmic ray interactions with the intracluster medium and the 
annihilation of dark matter. We report on a search of measurements from the EGRET 
instrument aboard the Compton Observatory for gamma-ray emission from clusters of 
galaxies. This cluster survey includes a preliminary analysis of data from both pointed 
observations and the EGRET all-skv survey. Current upper limits are given and the 
cosmological implications of the observations are briefly discussed. 

INTRODUCTION 

The diffuse gamma-ray emission in our galaxy is thought to be due to interactions 
between cosmic rays and the interstellar medium and baa been used to probe the 
properties of cosmic rays in the galaxy (Bloemen 1989; Fichtel c i at, 1989, and references 
therein). The presence of a relatively dense (n^ I0“ 3 cm" 3 ) intracluster medium (1CM) 
in many clusters of galaxies means that similar mechanisms can generate gamma-rays 
at those extragalactic sites. Many rich clusters of galaxies also contain active galaxies 
that are known to be powerful sources of relativistic electrons. For these, and other 
reasons that we discuss below, EGRET observations of clusters of galaxies can place 
interesting constraints on high energy phenomena at extragalactic sites. 

GAMMA RAYS FROM CLUSTERS OF GALAXIES 


Estimating the density of cosmic rays 

In models where the cosmic rays in the intracluster medium (1CM) are generated in 
the cluster galaxies, the number in a cluster it the number leaked from the constituent 
galaxies. If the generation of cosmic rays is proportional to optical luminosity, and we 
assume our Galaxy is typical, then 


A'c / = A, a , 


tC 


0 ) 


where t H is the Hubble time, t c is the confinement time of C 06 mic rays in the galaxy, A f cl 
is the number of cosmic rays in the cluster, and N ga i * 10 s * (Stecker and Tylka 1989), 
the number currently in our Galaxy. The L’b are the respective optical luminosities. 

Since the same emission mechanisms will be available in the cluster bj in the galaxy, 
we can calculate the gamma-ray luminosity of the duster as 
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where we have assumed that the cosmic rays and gas are distributed in a King model 
of index » 3 in the ICM. 

For a cluster at 100 Mpc this is still a factor of about 6 too weak to be detected in 
a 500,000 second observation if the ratios are unity. However, given the considerable 
uncertainties in many of the quantities (particularly in the confinement time) that go 
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into this luminosity, gamma rays from leaked cosmic rays may still be delectable by 
EGRET, particularly if the Hubble constant is large (so that clusters are relatively 
nearby). Even if individual clusters are not visible, a sample of clusters may show & 
significant excess of positive deviations. 

Cluster Rad\o Hales 


Radio observations indicate that relativistic electrons are present in the ICM of at 
least half a dozen clusters (Hanisch 1982). These dusters contain extended regions of 
low radio surface brightness that are associated with the cluster as & whole rather than 
individual active galaxies. A well studied example is in the Coma cluster where the 
halo has a scale size of roughly ] Mpc and a steep power-law spectrum (a ^ 1.3, e.$ M 
Kim et al. 1990). The radio measurements are best explained by synchrotron emission 
from relativistic electrons in an intracluster magnetic field. 

The spatial extent of Coma-type halos has led to difficulties for understanding their 
origin. It is natural to assume that the halo is powered by electrons that leak out 
of active radio galaxies. However, Jaffe (1977) has shown that these leakage models 
have difficulty accounting for the size of the radio halos. The crux of his argument 
is that if the bulk drift velocity of electrons is limited to the Alfven speed (Wentzel 
1974) then the time for electrons to diffuse to the outer parts of the radio halo exceeds 
the synchrotron* Compton loss timescale by at least two orders of magnitude. So while 
some energetic electrons must certainly leak into the ICM from active galaxies, the bulk 
of the radiating electrons must be generated in some other manner. Two approaches 
to the generation of relativistic electron* within the cluster ICM have been explored. 
The first approach, employed by ‘‘primary electron’’ models, assumes that the electrons 
are accelerated in situ by either shocks or turbulence in the ICM (Schlickeiser et al. 
1987). The second approach, employed by “secondary electron* models, assumes that 
the radiating electrons are generated when leakage cosmic ray protons, whose much 
longer radiative lifetime* allow them to fill the cluster volume, collide with thermal 
ICM nuclei (Dennison 1980). i 

Inasmuch a* both types of modeli for Coma-type halos. require the presence of high- 
energy particles in the ICM, it is if interest to examine gamma-ray production by 
these particles as they interact with the thermal gas required by cluster x*ray obser- 
vations. By combining the ICM deniity measurements provided by -x-ray observation* 
and the shape of non-thermal electron spectrum provided by the radio measurements, 
one can calculate the gamma-ray flujx given an assumed magnetic field strength. For 
the Coma cluster, the absence of a hard component of Compton x-rays limits the 
minimum field strength to 0.1 ^iG (Rephaeli et at. 1987). For primary electrons mod- 
els this yields a predicted maximum nonthermal bremsstrahlung flux at 100 MeV of 
~ 1 0~ e cm“ 2 .s" 3 M eV~ J . which would be difficult to detect with EGRET. On the oilier 
hand, secondary electron models predict a detectable pion decay flux from Coma if t)u* 
field strength is less than O.SpGauss (K/estrand 1990). A detection of Coma by EGRET 
would therefore be strong evidence far a reservoir of cosmic-ray protons in the ICM. 

Mtiagt dactic cosmic rays 

actic origin, then the flux of cosmic rays in dusters 
nple geometric arguments 


If local cosmic rays are of metagal 
should be easily detectable. Using si 
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where the densities are of the ICM versus ISM and s is the scale of the system. Sino* 
the product of these ratios is about unity, and the clusters fill a significant fraction <* 
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an EGRET resolution element, clusters should be visible. 

Annihilation 0} dark matter 

of dart* mr t e ( rally i! )elieve<1 th *- mc ? 1 of the matter in clusters is contained in some form 
S l ? r Wl ?°! e n ‘ tUfe '* unknow " . There are many possibilities for the nature 
gamma ravs^rd^otb & °> exot ' c P«-rticlee which may occasionally decay to 

?QR7 B»l y ! j c h * p ? du< ; t8 durm « enc °unters (see for example Silk and Bloemen 

Z o^trn iMv C f ^ il h Iiktly that loCal G a ' a ^ constraint. onT 

chance th^u/h d ' out a s 'S nlficant contribution from this process, there is some 
cnance that such decays could be measurable in clusters. 

cluster. W are & g ®° metric fact0r the reIa ‘ ive «urface brightnesses of the Galaxy and 
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where the n’s are the densities of the exotic particles and the s’s are the scale lengths 

distribution of f ac ^ or< ° f ° rde , r Unity are indetcrmilu ‘« without a knowledge of the 
maUCr “ d C ° Uld C ° nCeiv4bly fftVor detecli0 " cluster,® 
the 4,i> °r " P&rt,C,pate in the formati ° n of galaxies to the extent that 
> ^ iMpe — ‘ ban the cluster emission of annihilation products would 

^enhanced 1 f^tn '^cSVnn V S ' U th ' av L er& «®,^ uara of the density of exotic particles 
dusters lh d en r f 00 0r , m0f t ln ! he Ga!ayy above the corresponding average in 
toune extend - " t ’ U ' rf#e * br ]S h ' n f 88 wi » be g™ a ‘er than clusters (and will of 

matter Still riven t\ 7^* ^, ucb arger enhancements are seen in the visible 

its datftrtmn * T ^ 18 un ^ely to participate in dissipative proceasei 

“^;s'3r b ' r” p „s. mr “ ,ible * nd wi,b - ,iu " — — »*■-; 

THE OBSERVING PROGRAM 

archive ^ r"*' 3 "? , of both archival information from the CGRO 

of gamml ravs ^ The nK?" ° f ^ m08t Promising candidates for the detection 

ot gamma rays. The observing program was begun in Phase .2 with severs! nk 

Tmuch dM n ° d8 h 201 &nd 2 ° 2 < ^ r ‘ her observations are scheduled for Phase 3 including 
a much deeper observation of the Coma cluster. 8 

RESULTS 

in TnbtT' ! A?irj, ? '““'"T.? , " hniq ”' h “ b “” ‘"' <1 10 O' 1 '™ 1 "' th ' »«* Ita>>« 
n »w, S i , nd the 80urce re ? ion is chosen for background. The counts 

‘ daCkgr ° U regl ° n multiplied by the ™lio of the total exposure in the source 

the source ^rion^An J" * b * background a " nulus and subtracted from the counts in 
t e source region. Any excess counts are then attributed to the source. The source 

aridTn 'outer 7 H ° ? be back * round annulus has an inner radius of 9,5° 

used in tk f | d ' U3 °l L 3 fl ' 0n,y hlgb quality photons with energy > 70 MeV were 
region Note^at'thJ.v flUX ' 8 computed uain « lhe av erage exposure over the source 

by < E >- ( r p/o wn Wo ag t EXP T Ure in the region “ f8lited ‘he total exposure 
i/thefnf 7 whete ‘he integral over the solid angle of the source region 

is the total exposure in that region. 6 

with^the AGN X 4C?fi?l ° f , A2 ( 199 ' none of the ^ters was detected. A2199 is confused 
! " about a de « ree awa y ( Ma “° x */. 1992). Since the emission 

attributed . T 8tr ° n 2, vanallon8 a timescale most of the fiux must be 
attributed to a non-extended source. The fiux we obtain for this source agrees well 
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Cluster 

D 

(Ho =* 50) 
Mpc 

Viewing 

Periods 

Flux Limit 
(>100 M«V) 
10- T (cur l cm' , ) 

Luminosity Limit j 
(>100 MeV) j 
(10 44 erga *“ l ) 

A400 

139 

21,0 

1.9 

1 

4.0 

A426(Peneue) _ 

110 

15.0 



A496 

192 

29.0 1 

5.0 

20. 

A 1060 

68 

0.7 

4.5 

2.3 

A1S67 

“T5T 

3.0,4.0,11.5 


20. 

A1656(Coma) 

140 

3.0,4.0,11.0 

1.3 

2.8 

A2&22 

330 

9.5,201.0,202.0 | 

2.3 

- . , — — — 

27. 

A2Q52 

210 

24.0,24.5,25.0 

2.0 

9.6 

A2063 

200 

24.0,24.5.25.0 

’ To 

17. 

— A~' ~ 

A2162 

19Cp 

9.5,201.0,202.0 

O 71 

24. 

A2199 

180 

9.5,201.0,202.0 

HP 

35. 

- 1 7ft 

A2255 

480" 

9.5,201. 0,202"0 

2.8 

70. 

A 2256 

360 

‘9.5,201.0,202.0 

1.5 

21 , 

A2319 

336 

2.0, 7.0 

Ll 


A 24 62 

420 

9.0,13.5,19.0 

1.5 

29. 

A2806 

r§o" 

' 9.0,10.0,13.5 

1.4 

3.9 

A2870/A2877 

150 

9.0,10.0,13.5 

6.2 

15. 

A3526 

66 

12.0,14.0 

^ 3.2 

1.5 

A3565 

66 

12.0 

^ 2.0 

0.94 

A3574 

84 

12.0 

5.4 

4.1 

A4038 

170 

9.0,13.5 

1.4 

4.4 

A4049 

u__ ^ 

9.0,13,6 

1.7 

5.3 


' Confuted with Galactic flux. 


5 Confuted with 4C+38.41. 

Luminosity limits ataume a photon power lew index of 2.2. AU limits are 2o fluctuation. 
..Table 1. Clusters observed. 
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with t ha-t of the EGRET team. Upper limits on the flux of clusters which we have 
analyzed are described in Table 1. The simple aperture photometry approach we have 
used fails for A 426 and A2319 which are at low galactic latitude. The background for 
these sources must be handled in a more sophisticated fashion which will be possible 
with the likelihood technique. 

DISCUSSION 

The investigations we have so far presented are only preliminary. We anticipate 
that limits perhaps an order of magnitude deeper can be made with better analysis. 
Sophisticated likelihood analysis programs should soon become available and deeper 
observations of Coma will be made in Phase 3. Nevertheless, some of our results are 
already interesting. 

The non-detection of clusters drives another nail in the coffin of metagalactic origin 
for cosmic mys and supports the findings of Sreekumar tt ai (1902) who looked at 
emission from the Magellenic clouds. 

The limits on the emissivity of the Coma cluster increase the lower limit of the 
intracluster field strength by a factor of five to — Q.5pG if the radio halo is powered 
by secondary electrons. Using this lower limit in a simple turbulent cell model for 
the cluster field structure, the observations of the Faraday depth for polarized radio 
sources (Kim ct ai 1990) would require at least 4 turbulent cells per gas core radius. 
Wc estimate that the deeper exposure for Coma that should occur during the Phase 
3 observations will allow us to detect pion gamma-rays even if the field strength is 
lpG. Failure to detect Coma at that sensitivity level would force secondary electron 
models into a rather small region of parameter space. 

Using the simple model discussed above the confinement times for cosmic rays in the 
cluster galaxies must be > 10 T years. This is consistent with current estimates of the 
confinement times of several times 10 T yearB from analysis of the isotopic composition 
of local cosmic rays (e.^., Webber tt ai 1990). It should be possible to get very close 
to those limits with longer observations and enhanced analysis techniques. 
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